Rupture Kinematics of the 2005 M_w 8.6 Nias–Simeulue Earthquake from the Joint Inversion of Seismic and Geodetic Data

The 2005 M_w 8.6 Nias–Simeulue earthquake was caused by rupture of a portion of the Sunda megathrust offshore northern Sumatra. This event occurred within an array of continuous Global Positioning System (GPS) stations and produced measurable vertical displacement of the fringing coral reefs above the fault rupture. Thus, this earthquake provides a unique opportunity to assess the source characteristics of a megathrust event from the joint analysis of seismic data and near-field static co-seismic displacements. Based on the excitation of the normal mode data and geodetic data we put relatively tight constraints on the seismic moment and the fault dip, where the dip is determined to be 8° to 10° with corresponding moments of 1.24 x 10^(22) to 1.00 x 10^(22) N m, respectively. The geodetic constraints on slip distribution help to eliminate the trade-off between rupture velocity and slip kinematics. Source models obtained from the inversion of various combinations of the teleseismic body waves and geodetic data are evaluated by comparing predicted and observed long-period seismic waveforms (100–500 sec). Our results indicate a relatively slow average rupture velocity of 1.5 to 2.5 km/sec and long average rise time of up to 20 sec. The earthquake nucleated between two separate slip patches, one beneath Nias and the other beneath Simeulue Island. The gap between the two patches and the hypocentral location appears to be coincident with a local geological disruption of the forearc. Coseismic slip clearly tapers to zero before it reaches the trench probably because the rupture propagation was inhibited when it reached the accretionary prism. Using the models from joint inversions, we estimate the peak ground velocity on Nias Island to be about 30 cm/sec, an order of magnitude slower than for thrust events in continental areas. This study emphasizes the importance of utilizing multiple datasets in imaging seismic ruptures

Rupture Process of the 2004 Sumatra-Andaman Earthquake

The 26 December 2004 Sumatra-Andaman earthquake initiated slowly, with small slip and a slow rupture speed for the first 40 to 60 seconds. Then the rupture expanded at a speed of about 2.5 kilometers per second toward the north northwest, extending 1200 to 1300 kilometers along the Andaman trough. Peak displacements reached ~15 meters along a 600-kilometer segment of the plate boundary offshore of northwestern Sumatra and the southern Nicobar islands. Slip was less in the northern 400 to 500 kilometers of the aftershock zone, and at least some slip in that region may have occurred on a time scale beyond the seismic band

Characterizing fluid flow paths in the Hellisheidi geothermal field using detailed fault mapping and stress-dependent permeability

The Husmuli zone of the SW-Iceland Hellisheidi geothermal field is currently being used for re-injection of geothermal fluids and geothermal CO2 for its permanent storage in the form of carbonate minerals. A fully coupled hydro-thermo-mechanical numerical model was employed to investigate the coupled impacts of these complex processes on the calibration of fluid flow paths, which can have significant implications for the long-term performance of this subsurface reservoir. Employing a combination of high-resolution fault mapping with laboratory measurements of stress dependent permeability coupled into a dual porosity field-scale model, the flow paths were calibrated using results of tracer tests performed at the site using stress-dependent permeability tensors. Although vertically extended faults are the primary fluid flow paths, fractures connecting the faults can play an important role in fluid transport. As the upward flow streamlines manifest, deep geological layers can also deviate the fluid flow towards the shallower layers provoking the vertical flow of geothermal fluids. This highlights the sweet spot for sustainable flow and heat extraction in vicinity of faults intercepting the geological layers at depth of 1100 m. It was also shown that the inclusion of the geomechanical calculations in the history matching of the tracer test could lead to changes in arrival time and peak of the tracer profiles. Results of an independent tracer test were used to validate the model and to demonstrate the predictive capability of the calibrated model. This verifies the consistency of our methodology to incorporate the stress-dependent permeability. The results of this comprehensive modelling study provide insight into the likely fluid flow paths, which can have profound impact on the evaluation of various processes such as CO2 mineralisation taking place in Hellisheidi geothermal reservoir

Induced seismicity risk analysis of the hydraulic stimulation of a geothermal well on Geldinganes, Iceland

The rapid increase in energy demand in the city of Reykjavik has posed the need for an additional supply of deep geothermal energy. The deep-hydraulic (re-)stimulation of well RV-43 on the peninsula of Geldinganes (north of Reykjavik) is an essential component of the plan implemented by Reykjavik Energy to meet this energy target. Hydraulic stimulation is often associated with fluid-induced seismicity, most of which is not felt on the surface but which, in rare cases, can be a nuisance to the population and even damage the nearby building stock. This study presents a first-of-its-kind pre-drilling probabilistic induced seismic hazard and risk analysis for the site of interest. Specifically, we provide probabilistic estimates of peak ground acceleration, European microseismicity intensity, probability of light damage (damage risk), and individual risk. The results of the risk assessment indicate that the individual risk within a radius of 2km around the injection point is below 0.1 micromorts, and damage risk is below 10-2 for the total duration of the project. However, these results are affected by several orders of magnitude of variability due to the deep uncertainties present at all levels of the analysis, indicating a critical need in updating this risk assessment with in situ data collected during the stimulation. Therefore, it is important to stress that this a priori study represents a baseline model and starting point to be updated and refined after the start of the project

Rupture process of the 2004 Sumatra-Andaman earthquake.

The 26 December 2004 Sumatra-Andaman earthquake initiated slowly, with small slip and a slow rupture speed for the first 40 to 60 seconds. Then the rupture expanded at a speed of about 2.5 kilometers per second toward the north northwest, extending 1200 to 1300 kilometers along the Andaman trough. Peak displacements reached approximately 15 meters along a 600-kilometer segment of the plate boundary offshore of northwestern Sumatra and the southern Nicobar islands. Slip was less in the northern 400 to 500 kilometers of the aftershock zone, and at least some slip in that region may have occurred on a time scale beyond the seismic band

Full-Waveform based methods for Microseismic Monitoring Operations: An Application to Natural and Induced Seismicity in the Hengill Geothermal Area, Iceland

Geothermal systems in the Hengill volcanic area, SW Iceland, started to be exploited for electrical power and heat production since the late 1960s. Today the two largest operating geothermal power plants are located at Nesjavellir and Hellisheiði. This area is a complex tectonic and geothermal site, located at the triple junction between the Reykjanes Peninsula (RP), the Western Volcanic Zone (WVZ), and the South Iceland Seismic Zone (SISZ). The region is seismically highly active with several thousand earthquakes located yearly. The origin of such earthquakes may be either natural or anthropogenic. The analysis of microseismicity can provide useful information on natural active processes in tectonic, geothermal and volcanic environments as well as on physical mechanisms governing induced events. Here, we investigate the microseismicity occurring in Hengill area, using a very dense broadband seismic monitoring network deployed in Hellisheiði since November 2018, and apply sophisticated full-waveform based method for detection and location. Improved locations and first characterization indicate that it is possible to identify different types of microseismic clusters, which are associated with either production/injection or the tectonic setting of the geothermal area